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DESCRIPTION 

Method of determining a structure of a moving object 

5 The present invention relates to the field of digital imaging. In particular, the present 
invention relates to a method of determining a structure of a moving object from an at 
least two dimensional data set, an image processing device and a computer program for 
an image processing device. 

10 The anatomy of the human coronary tree is widely consistent between normal subjects. 
It always consists of three main branches called left anterior descending (LAD), left 
circumflex (LCX) and right coronary artery (RCA). The main variation in topology to 
be found is a dominance of either left or right branch supplying the apex of the heart, as 
described in M D. Cerqueira et al. "Standardized Myocardial Segmentation and 

15 Nomenclature for Tomographic Imaging of the Heart" Circultation 105: 539-542, 2002 
which is hereby incorporated by reference, attaching more sub-branches to either side. 
Geometric properties of segments, however, may vary significantly between subjects, 
for example, due to differences in size and shape of the ventricles of the heart 

20 Document C. Chalopin, LB. Magnin, G. Finet "Automatic labeling of the coronary tree 
using a three-dimensional reference prior Model" Computers in Cardiology 25, p. 761- 
764.1998 discloses the generation of two dimensional reference models of either left or 
right coronary tree from a given three-dimensional model, each of which having a 
different projection angle. The best fitting reference model is used to identify a coronary 

25 main branch in geographic projections and to label the segments of the image. In other 
words, the best fitting reference model is used to identify the respective branches of the 
human coronary tree in an image. 

It is an object of the present invention to provide for an improved prediction of 
30 properties of a moving object such as, for example, the human heart 
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The above object may be solved by a method of determining a structure of a moving 
object from an at least two-dimensional data set according to claim 1. According to this 
method, a model of the structure relating to the properties of interest is applied to the 
data set. Then, an adaptation of the model to the data set is performed Then, on the 
5 basis of the adapted model, a location of at least one portion of the structure of the 
moving object is estimated 

In other words, with respect to the example of the coronary tree of a human heart, a 
model of the coronary tree, which is, for example, a statistic deformable model based on 

10 preceding measurements, is applied to, for example, a two-dimensional image, such as 
an x-ray angiography image. Then, this model, which may be a three-dimensional 
model, is adapted to the two-dimensional image representing the data set. This 
adaptation may, for example, be performed by a parameter variation. Then, a location of 
at least one portion of the structure, for example of a portion of the coronary artery with 

15 a total occlusion, which renders the vessel invisible, since no contrast agent is supplied 
to this region, is predicted by using the adapted model. 

Advantageously, this exemplary embodiment of the present invention may allow for 

improveddia^osticso^ fbrexam 

20 detection of side branches of the coronary tree, since, on the basis of the model, a 

location of at least one portion of the coronary tree or the artery tree may be predicted 

According to another exemplary embodiment of the present invention as set forth in 
claim 2, the adapted model is overlaid onto the image of the coronary tree such that an 

25 imsge may be displayed showing the actual measured image and the overlaid model,. 
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According to another exemplary embodiment of the present invention as set forth in 
claim 3, parameters of the model are adapted on the basis of a similarity of the model to 
the structure, which may allow for a simple and efficient adaptation of the model. 

5 According to another exemplary embodiment of the present invention as set forth in 
claim 4, an image quality of an image displayed may be improved by superimposing 
images of the object of interest taken at different points in time or with different 
projections. For superimposing images taken at different points in time or projections, 
the adapted model may be used as a reference, allowing for a registration of these 
10 images. 

According to another exemplary embodiment of the present invention as set forth in 
claim 5, the model is a deformable model and the adaptation of the model is performed 
by an energy minimization of an internal and an external energy of the model, which 
1 5 allows for a fast adaptation of the model. 

According to another exemplary embodiment of the present invention as set forth in 
claim 6, the model is a statistical model of a coronary tree of a human heart and the data 
set relates to x-ray angiography data. 

20 

According to another exemplary embodiment of the present invention as set forth in 
claim 7, an image processing device is provided, allowing to estimate a location of at 
least one portion of a structure of an object of interest, which, for example, is not 
discernible in an image representing a two-dimensional data set. This estimation is 
25 performed on the basis of a model of the structure of interest adapted to the image, i.e. 
the data set 

Further exemplary embodiments of the image processing device are provided in claims 
8and9. 

30 



-4- 



PHDE030402 EPP 



According to another exemplary embodiment of the present invention as set forth in 
claim 9, a computer program is provided allowing for an estimation of a location of at 
least one portion of a structure of a moving object of interest by using a model of the 
structure of interest adapted to the measured data set The computer program may be 
5 written in any suitable programming language, such as C++ and may be stored on a 
computer readable medium, such as a CD-ROM. However, the computer program 
according to the present invention may also be presented over a network such as the 
WorldWideWeb, from which it may be downloaded into the working memory of a 
processor. 

10 

It may be seen as the gist of an exemplary embodiment of the present invention that a 
location of a portion of a structure on or of a moving object, which, for example, is not 
discernible in an image representing the original measured data set is predicted by 
using a model adaptation. For this, a model of the structure of interest is adapted to the 
1 5 structure of interest in the original data set Once the model is adapted to the image, 
portions (a location, a dimension or a direction thereof) of the structure not discernible 
or visible in the original measured image may be indicated by the model which may be 
displayed in an overlaid fashion on, for example, a display. 

20 These and other aspects of the present invention will become apparent from and 
elucidated with reference to the embodiments described hereinafter. 

Exemplary embodiments of the present invention will be described in the following, 
with reference to the following drawings: 



25 
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Fig. 3 shows a schematic representation of a coronary model according to an exemplary 
embodiment of the present invention. 

Fig. 4 shows a projection of a geometric coronary mean model with branches 
corresponding to Fig. 3 according to an exemplary embodiment of the present 
invention. 

Fig. 5 shows a sub-set of branches overlaid on an image and a corresponding area at a 
well-fitted model, according to an exemplary embodiment of the present invention. 

Fig. 6 shows a flow-chart of another exemplary embodiment of a method according to 
the present invention. 

Fig. 7 shows images of a coronary vessel tree generated in accordance with the present 
invention. 

Fig. 8 shows a flow-chart of another exemplary embodiment of a method according to 
the present invention. 

Fig. 1 shows a simplified schematic representation of an exemplary embodiment of an 
image processing device in accordance with the present invention. In Fig. 1, there is 
shown a central processing unit (CPU) or image processor 1 for applying a model of the 
structure of interest to the original measured data set, for performing an adaptation of 
the model to the data set and for performing an estimation of a location of at least one 
portion of the structure by using the adapted model. The image processor 1 is connected 
to a memory 1 for storing the data set, for example, a plurality of images relating to x- 
ray angiography images. The image processor 1 may be connected by a bus-system 3 to 
a plurality of periphery devices or input/output devices which are not depicted in Fig. L 
For example, the image processor 1 may be connected to an x-ray scanner. However, 
the image processor 1 may also be connected to an MR device, to a CT device, to an 
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ultrasonic scanner, to a plotter or a printer or the like via the bus-system 3. Furthermore, 
the image processor 1 is connected to a display such as a computer screen 4, for 
outputting information and/or images to a user. Furthermore, a keyboard 5 is provided, 
connected to the image processor 1, by which a user or operator may interact with the 
5 image processor 1 or may input data necessary or desired for the segmentation process. 

Fig. 2 shows a flow-chart of an exemplary embodiment of method according to the 
present invention, which may be used for operating the image processing device 
depicted in Fig. 1. 

10 

In the following, the present invention will be described with respect to the 
determination of the anatomy of the human coronary artery tree. In this connection, the 
present invention is particularly suited to improving and speeding up the diagnostics of 
coronary artery defects and their treatment Thus, for example, the present invention 
1 5 may be applied as an intervention support system for 2D-3D x-ray angiography. Also, 
the present invention may be applied in conjunction with catheter intervention. 

As may be taken from Fig. 2, in step SI, an image is acquired, for example, by x-ray 
^ - angiography: This image may relate to a-particular phase of-the heart and to a-particular. . 

20 projection angje. Then, in step S2, the model, which may for example be a statistic 

model projected into the plane of the image, is used to measure a distance between parts 
of the model projected into the plane of the image and corresponding structures in the 
image. This distance measured in step S2 is used to determine a similarity value, in step 
S3, reflecting a similarity of the coronary artery tree depicted in the image and the 

25* model to be adapted to the actually measured coronary artery tree in the image. In case 
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indicated that there should be an artery portion, but in the actually measured image there 
is no artery portion, the user may assume that, for example, there may be an occlusion 
in the artery, such that no contrast fluid in the blood makes the artery visible in this 
region. 

5 

As may be taken from steps S4 and S5, parameters of the model are varied. Then, as 
may in particular be taken from S5, positions and connections of the model may be 
varied and parameters such as heart phase q>, the projection angle p of the model and an 
individual derivation from the mean model given by <p and p may be varied, which is 
1 0 indicated in step S6, allows a prediction of model positions by calculating a projection 
of model positions. In other words, in step S6, the three-dimensional model varied and 
adapted in the preceding steps, is projected into the image plane of the image acquired 
in step SI . Then, as indicated in step S7, the projected image is used to re-enter the 
iteration in step S2. 

15 

The coronary model used according to an exemplary embodiment of the present 
invention may, for example, be taken from Figs. 3 and 4. Fig. 3 shows a schematic 
coronary model (topology only) with LAD, LCX and RCA as well as some sub- 
segments according to an exemplary embodiment of the present invention. Fig. 4 shows 
20 a projection of the geometric coronary mean model with branches corresponding to Fig. 
1 according to an exemplary embodiment of the present invention, as it may be 
achieved in step S7. 

The geometry of each branch of the geometry model is represented parametrically in 
25 such way that it expresses the position of a branch both in global and local coordinates. 
In the local coordinates, the position of a branch is defined relative to another branch, 
for example, to a parent branch. The parameters of the model are modeled statistically 
such that a given parameter setting is assigned a probability value. A specific potential 
geometry of the coronary tree, as indicated in steps S2 to S7 in Fig. 2 may be predicted 
30 by a parameter setting. Movement and deformation of the branches caused by heartbeat 
and by respiration are covered by the model. For the heartbeat, this means that geometry 
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predictions are possible for any heart phase. To illustrate the use of the coronary model 
in more practical terms, an exemplary embodiment of a parameterization and a 
sampling is described in the following. 

5 According to this exemplary embodiment of the present invention, each branch is 
assigned a set of sample positions s bJ . Each branch b is interpolated through its 

sampled position by a spline. The model holds a statistical description of reached 
sample position, for example, a mean value s bi and covariance matrix. These values 
may be referenced to a global coordinate system, or they may be positioned relative to 

10 the corresponding parent branch (i.e. global or local coordinates). Also, information 
about lumen diameter of each branch may be provided by the model by l u . The 
deformation of the coronary tree is covered by a variation of sample positions over the 
heart cycle <p = 0..2rc. Similar to the spatial sampling by sample positions, some 
temporal samples s^u with t=L.m are part of the model, where a temporal interpolation 

15 of the distribution of positions is provided between the m temporal samples. Thus, the 
model according to this exemplary embodiment of the present invention may be used to 
predict the position of each branch b at a given heart phase <p by a statistical description 

--of the sample positions be "described ~ " 

in further detail in the following. 



For a given heart phase q> g and a given projection plane orientation p gt the projection of 



20 
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However, individual properties of the patient's heart may deviate from the mean model 
s bJ j g . Also the given parameters like the projection orientation p g or especially the heart 

phase (p g may not precisely reflect the image acquisition situation. Thus, the model 
5 should be allowed to warp itself in order to find a good congruence between model 
branches and visible branch position in the angiographic image. 

According to an exemplary embodiment of the present invention, this may be done by 
applying a deformable model approach where the model parameters s 6 ^ m , q> m , and p m 

10 are varied to optimize a similarity measure a = d m + f m (l) (see steps S2 S3 in Fig. 2). 

This similarity measure may be given by an image feature term f m (l) which determines 
how good the projected model branch positions agree with features in the image that 
correspond to coronaries (see step S2 in Fig. 2). A good proposal for f m (l) is a measure 
15 of the directed gradient at the predicted vessel boundaries. It gives a tnayimiim value, 
when predicted coronary positions in the projection coincide with real vessel walls in 
the image of the original data set. At both sides of the spline through all sample 
projections ofa given branch P s [$ g *P g )=\ s b g M 0 M p g \ and with a distance /^/ to it, a 

strong image gradient orthogonal to this line is expected. It is measured at a couple of 
20 sample points and all values may be averaged. 

According to another exemplary embodiment of the present invention, f m (/) 
corresponds to the sum of gray values within the area of all predicted projected 
coronaries, assuming that coronaries appear bright in the image. For each position in 
25 P s , p g ) = { • M P9 j given by the actual parameter setting of the model, a circular 

area Cbj with diameter l bJ is taken. The unification of all areas C/may be given by 
U = ^C bl . The average image intensity on this selected part of the image yields the 

feature term 
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J/(x) dx 

f =M 

W Jldr ' 

u 

Fig. 5 shows a sub-set of overlaid on the image and the corresponding area U (union 
of circular areas) at a well-fitted model in accordance with an exemplary embodiment of 
the present invention. As may be taken from Fig. 5, Fig. 5 shows an example of the area 
5 U where the image intensity is measured By using, for example, steerable filters on the 
image at the modeled position, both image gray value properties and image gradient 
properties may be combined. However, as a downside, using such steerable filters 
increases computational efforts required. 

10 The distance measure d m of the parameter setting to the modeled distribution is 
considered in the similarity measure. There must be distance measures for each 
parameter or parameter set: 

1 5 for each b and each z, 

hJa \=\U. —A II 

_.. -<p\rm/ [prm-- rgp 

and 

<t p (P m ) = \\Pm-Pg\\ 

where [|x|| is a suitable distance measure in the domain of the parameter. As the model 
20 holds a distribution d s ^ for each parameter s b i ^ also a statistical distance measure 

fe?>5"» TVi3 u3ii SuulGuj £T dlST37TC£ iJiSTSLiSO CGnSItferS "ui3 CGVST1 3HC2" jjjSjjl ju Ji IS SurOiICSbls. A 
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This term introduces constraints by the a priori knowledge about the expected shape of 
the vessel tree projection and will penalize parameter configurations that might fit well 
with the data but that are very unlikely to reflect an existing imaging situation. 

5 According to another exemplary embodiment of the method according to the present 
invention, a segmentation step to find the main branches in the angiographic image 
precedes the adaptation. It yields a measure of the center lines of successfully 
segmented coronaries called S. The subsequent adaptation is computationally less 
extensive, because no feature term has to be considered during the adaptation and no 
1 0 access to the image is required anymore. Instead of the image feature term f m (/) a 
distance measure £ m {S) on the segmentation result is used: j = cf m + <5 m (S) . S m (S) 
determines the distance between s bJA _ • M D and the center line of the segmented 

coronaries in the image. This variant also allows to use different model candidates and 
to decide, which one fits best. There could be for instance a model for left-dominant 
1 5 type and one for the right-dominant type. 

Instead of a global measure for the whole coronary tree as those introduced above, also 
a part of the model representing a single branches or sub-trees may be adapted to an 
image. The position of such a sub-tree may be modeled statistically with respect to its 
20 parent branch. Thus, the given or previously estimated position of this parent branch 
may be used to predict the position of its desired sub-tree in the image. 

Instead of using a set of template models, each with a given geometry, a deformable 
model is used that allows its geometrical parameters to be varied within a certain range 
25 of values. Once it is adapted to fit the available data, a geometric prediction 3D and 
even in 4D (space + time) is available for the whole coronary tree. The range of 
parameters is determined by the modeled statistical distribution of the parameters, but 
also by the result of the feature term f m (l) applied to the image - or S m (S) applied to 
the extracted center lines respectively. 
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Thus, the model may be adapted to an individual image wherever evidence is 
sufficiently available from the image. In image regions where evidence given by f m (l) 
is low, a priori information d m from the model is dominating the adaptation result This 
5 advantageous property of the model and the application of the model in accordance with 
the present invention allows a well-founded estimation of branch position taking into 
account available knowledge about proximate branches. In the same way, a trajectory of 
a branch in time may be predicted from one or more images taken at different heart 
phases. The geometry at an unknown heart-phase may be estimated from the model 
10 taking into account available data from those heart phases for which images are 

available. Thus, advantageously, for example, patients with chronic total occlusions 
(CTO) in coronary arteries, the whole vessel branch beyond the CTO is usually not 
visible in angiographic images, because no contrast agent is transported there, can be 
diagnosed more accurately and easily. 

15 

Fig. 6 shows a flow-chart of another exemplary embodiment of a method according to 
the present invention. This exemplary embodiment of the present invention is in 
particular suited, for example, for diagnosis of CTOs in coronary arteries where the 

20 because no contrast agent is transported there. 

The two-dimensional images of angiographic images makes it hard to the cardiologist to 
predict the position and extent of coronaries. It is especially hard to differentiate an 
occluded ending (CTO) from an ending that extends further but orthogonal to the image 

0^ -a%1 crpj a Tr ca pr vniTi «* QT^ ra r>r^jjyn ■gj -ieryn nf ggrygrgl vvroj ra^g TTVvrn A? i^u i>rp?if > crr^nrl ca<s> ja ii g? tolle r 
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a branch indicated as missing. This will be described in further detail with reference to 
Figs. 6 and 7. 

Steps SI to S7 of Fig. 6 correspond to steps SI to S7 in Fig. 2 and thus, to avoid 
5 unnecessary repetition, reference is made to Fig. 2 for the description of steps SI to S7. 

In contrast to the method of Fig. 2, from step S7, the formally three-dimensional model, 
which has been adapted to the image acquired in step SI, i.e. which has been adapted to 
feature in the images acquired in step SI and which has been projected into the image 

10 plane of the image acquired at step SI, is overlaid in step S10 onto the image acquired 
in step SI (an interruption of the adaptation circle of steps S2 to S7 to branch off to step 
S10 may be decided in step S3 when a sufficient similarity has been reached). Then, the 
image is displaced in step SI 1 to an observer. Due to this, the image visualized in step 
S10 shows the coronary tree or the respective arteries discernible, i.e. visible in the 

1 5 image and the overlaid adapted and projected model. Due to this, for example, a CTO 
case becomes obvious by a branch indicated by the model, but missing in the image. 

This is depicted in Fig. 7. Fig. 7 shows three different images generated in accordance 
with an exemplary embodiment of the present invention. The first image indicated in 

20 Fig. 7 by a shows a main artery (black), which is clearly visible in the image. However, 
as may be taken from image a, a proximal part of a coronary vessel tree is visible in this 
angiographic image (black) and a distal part, which is not visible in the angiographic 
image and which lies beyond a total occlusion is predicted by the model (white). As 
may be taken from image b in Fig. 7, the main branch has an occlusion and would 

25 usually not be visible in the angiographic image. However, the occluded part of the 
main branch is indicated by the model. Image c shows a functional occlusion. 

In other words, for CTO diagnostic imaging, the model is adapted to an angiographic 
image in the way explained above. A simulated projection of the model is then overlaid 
30 with the angiographic image. A fused visualization of a measure (the image) and the 
expected normal finding (the warped model) immediately allows the cardiologist to 
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diagnose a CTO (see Fig. 7) at positions where a predicted branch is missing in the 
image. Especially in screening situations and in cases without previous indication of 
CTO, it appears hard to visually scan all relevant positions. By using the present 
invention as indicated above, finding CTOs may be significantly facilitated because 
differences catch the eye. 

The above explained exemplary embodiment of the present invention may, for example, 
also be applied in conjunction with catheter intervention, where a catheter is to be 
navigated through the vessel to an occlusion. In particular, in cases such as depicted in 
image a (Fig. 7), it is usually hard to find the entry point of a side branch having an 
occlusion, since the side branch itself is not discernible in the measured image. The 
visualization of the predicted position, as shown in image a (Fig. 7), allows an improved 
navigation of the catheter. 

Fig. 8 shows a flow-chart of another exemplary embodiment of a method according to 
the present invention. As steps SI to S7 of the method depicted in Fig. 8 correspond to 
steps SI to S7 of the method depicted in Fig. 2, reference is made to Fig. 2 for a detailed 
description with respect to steps SI to S7. 

After a sufficient similarity has been reached, (which may be determined in step S3) the 
method may continue from step S5 to step S21. In step S21, branching off from step s5, 
a prediction is made by calculating projections of model positions for other heart phases 
<p and other projection angles p. Then, in the subsequent step S22, the warping 
parameters are determined and applied for overlaying images, which may allow for an 
image enhancement The operation performed in steps S21 and S22 will be described in 
farther detail in the following: 
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orientations are available and local image registration could be achieved The projection 
of a sample point of a given branch b g is given by • M fim where all model 

parameters are adapted to the image. These corresponding image positions are used to 
achieve a landmark-based registration between all the images. In order not to 
5 superimpose images from largely different viewing angles, only projection with low 
differences in p g - or p m respectively - should be chosen. The superposition of a line in 
the images that is expected from the adapted model to be the location of a branch makes 
the branch visible even though it was not clearly visible in the single images alone. This 
should work for smaller sub-branches that are hardly detectable and for the case of 
10 functional occlusion (Fig. 7) where only a certain amount of contrast agent is 

transported to. But this technique is even promising for totally occluded branches as it 
also enhances the poorly visible vessel wall. The true benefit from the model-based 
approach comes by aligning the model to the visible parts of the coronary tree and 
predicting the position of not sufficiently visible parts from that 

15 

According to exemplary embodiments of the present invention, a model of the coronary 
artery tree is adapted to the coronary artery tree visible in a measured image. The 
estimated position and orientation of cardiac structures reflected by the model may, 
according to an exemplary embodiment of the present invention, be used to calculate a 

20 set of parameters, for example, for the x-ray equipment. Such set of parameters may 
include position and orientation of x-ray generators and detectors in order to produce 
the cardiac standard projection image or in order to achieve an optimal x-ray set-up to 
image a given coronary artery. Furthermore, an iso-center for rotational cardiac x-ray 
imaging may be determined according to this exemplary embodiment of the present 

25 invention. By continuously feeding generated images during the examination to the 
model registration method, the accuracy of the adapted geometric model and the 
estimated position and orientation of the cardiac structures improve fiirther, especially if 
additional structures (for example, additional coronary arteries) are imaged. 
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Advantageously, this may allow for a very accurate setting of the x-ray equipment and a 
reduction in the x-ray dose applied to the object of interest, i.e. the patient This may 
also allow to reduce the amount of contrast medium applied to the patient. 

5 In spite of the fact that the present invention was described above with reference to the 
human heart and the determination of a coronary artery tree of the human heart, it is 
obvious to the skilled person that the above described present invention may also be 
applied to moving objects in general. In particular, the present invention may be applied 
to the determination of structures of a moving object from two-dimensional data sets, 
10 where, for example, the data set comprises images taken from different projections or 
taken at different points in time. 



15 
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CLAIMS 



10 



1 . Method of detennining a structure of a moving object from an at least two 
dimensional data set, the method comprising the steps of: 

applying a model of the structure to the data set; 

performing an adaptation of the model to the data set; 

estimating a location of at least one portion of the structure by using the adapted 

model. 



2. The method of claim 1 , 

wherein a first image is generated from the data set; and 
wherein, in at least the estimated location of the at least one portion of the 
structure, the adapted model is overlaid onto a second image which is based on the data 
15 set. 

3. The method of claim 1, 

wherein the data set comprises at least a third and a fourth image, wherein the 
third image relates to a first point in time or first projection and the fourth image relates 
20 to a second point in time or second projection, the first and second points in time and 
the first and second projections being different; 

wherein parameters of the model are adapted on the basis of a similarity of the 
model to the structure; and 

wherein the model is projected into a third projection of the first image. 

25 

4. The method of claim 3, 

wherein, for improving an image quality of a fifth image, the third image and the 
fourth image are superimposed; 

wherein, for superimposing the third and fourth images, the adapted model is 

30 used. 
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5. The method of claim 1, 

wherein the model is a deformable model; and 

wherein the adaptation of the model is performed by an energy minimization of 
an internal and an external energy of the model. 

5 

6. The method of claim 1, 

wherein the model is a statistical model of a coronary tree of a human heart and 
wherein the data set relates to x-ray angiography data. 

10 7. Image processing device, comprising: 

a memory for storing an at least two dimensional data set; and 
an image processor for determining a structure of a moving object from the at 
least two dimensional data set, wherein the image processor is adapted to perform the 
following operation: 
1 5 applying a model of the structure to the data set; 

performing an adaptation of the model to the data set; 

estimating a location of at least one portion of the structure by using the adapted 

model. 



20 8. The image processing device of claim 7, 

wherein a first image is generated from the data set; 

wherein, in at least the estimated location of the at least one portion of the 
structure, the adapted model is overlaid onto a second image which is based on the data 
set. 

25 
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wherein parameters of the model are adapted on the basis of a similarity of the 
model to the structure; and 

wherein the model is projected into a third projection of the first image, 
wherein, for improving an image quality of a fifth image, the third image and the 
5 fourth image are superimposed; 

wherein, for superimposing the third and fourth images, the adapted model is 

used. 

10. Computer program for an image processing device, wherein the computer program 
10 is adapted to cause a processor of the image processing device to perform the following 
operation when the computer program is executed on the processor: 
applying a model of the structure to the data set; 
performing an adaptation of the model to the data set; 

estimating a location of at least one portion of the structure by using the adapted 

15 model. 
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ABSTRACT 

In patients with chronic total occlusions (CTO) in coronary arteries, the whole vessel 
branch beyond the CTO is usually not visible in angiographic images because no 
5 contrast agent is transported there. According to an exemplary embodiment of the 

present invention, a model of a structure of interest of a moving object, i.e. for example, 
a model of the coronary artery tree of the human heart is applied to a data set or image 
of the object of interest. The model is adapted to features of the image. Then, a location 
of at least one portion of the structure is estimated by using the adapted model. 
10 According to an exemplary embodiment of the present invention, the adapted model is 
overlaid onto the measured image and displayed to a user, allowing, for example, for a 
determination of a CTO case, which becomes obvious by a vessel branch indicated as 
missing. 
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